MIRZAKHANI'S RECURSION RELATIONS, VIRASORO 
CONSTRAINTS AND THE KDV HIERARCHY 
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Abstract. We present in this paper a differential version of Mirzakhani's 
recursion relation for the Weil-Petersson volumes of the moduli spaces of bor- 
dered Riemann surfaces. We discover that the differential relation, which is 
equivalent to the original integral formula of Mirzakhani, is a Virasoro con- 
straint condition on a generating function for these volumes. We also show that 
the generating function for rjj and k\ intersections on M g , n is a 1-parameter 
solution to the KdV hierarchy. It recovers the Wittcn-Kontscvich generating 
function when the parameter is set to be 0. 



1. Introduction 

In her striking series of papers [18, 19], Mirzakhani obtained a beautiful recursion 
formula for the Weil-Petersson volume of the moduli spaces of bordered Riemann 
surfaces. Her recursion relation is an integral formula involving a kernel function 
that appears in the work of McShane [17] on hyperbolic geometry of surfaces. We 
have discovered that the differential version of the Mirzakhani recursion formula, 
which is equivalent to the original integral form, is indeed a Virasoro constraint 
condition imposed on a generating function of these volumes. 

Mirzakhani proves in [19] that her recursion relation reduces to the Virasoro con- 
straint condition as the length parameters of the boundary components of Riemann 
surfaces go to infinity, and moreover, it recovers the celebrated Wittcn-Kontsevich 
theorem of intersection numbers of tautological classes on the moduli spaces of sta- 
ble algebraic curves. Our result reveals that the Virasoro structure exists essentially 
in the Mirzakhani theory, and that it is not the consequence of the large boundary 
limit. 

The Virasoro constraint formulas for the generating functions of Gromov-Witten 
invariants of various target manifolds have been extensively studied in recent years 
[4, 8, 9, 22]. Although Mirzakhani's hyperbolic method does not immediately apply 
to these cases with higher dimensional target spaces, the Virasoro structure we 
identify in this paper strongly suggests that the Virasoro constraint conjecture of 
[5, 6] is a reflection of the combinatorial structure of building the domain Riemann 
surface from simpler objects such as pairs of pants or three punctured spheres. 

Although it is more than 15 years old, the Witten-Kontsevich theory [23, 15] 
has never lost its place as one of the most beautiful and prime theories in the 
study of algebraic curves and their moduli spaces. The theory provides a complete 
computational method for all intersection numbers of the tautological cotangent 
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classes (the ^-classes) defined on the moduli space M g , n of stable algebraic curves 
of genus g with n marked points. Recently several new proofs have appeared 
[21, 19, 14, 13]. We note that all these new proofs are based on very different ideas 
and techniques, including random matrix theory, random graphs, Hurwitz theory, 
representation theory of symmetric groups, symplectic geometry and hyperbolic 
geometry. 

The mystery of the Witten-Kontsevich theory has been the following question: 
where does the KdV equation, and also the Virasoro constraint condition, come 
from? Once we accept the Kontsevich matrix model expression of the generating 
function of all cotangent class intersections, then both the KdV and the Virasoro 
are an easy corollary of the analysis of matrix integrals. Thus the real question 
is: where do these structures appear in the geometry of moduli spaces of algebraic 
curves? 

The insight of some of the the new proofs [14, 13], which do not rely on matrix 
integrals but rather use the counting of ramified overings of P 1 , is that the KdV 
equation is a direct consequence of the cut and join mechanism of [10]. 

The proof [19] due to Mirzakhani utilizes hyperbolic geometry and has a markedly 
different nature from the others, whose origins are rooted in algebraic geometry. 
Mirzakhani's work concerns the Wcil-Petersson volume of the moduli space of bor- 
dered Riemann surfaces M. g , n (L), where L = (Li, . . . ,L n ) specifies the geodesic 
lengths of the boundaries of Riemann surfaces. Here the moduli space is equipped 
with the structure of a differentiable orbifold realized as the quotient of the Te- 
ichmiiller space by the action of a mapping class group. Mirzakhani shows that 
these volumes satisfy a recursion relation, and that in the limit L — > oo her recur- 
sion formula recovers the Virasoro constraint condition for the generating function 
of V'-class intersection numbers of M g ^ n , or equivalently, the generating function 
of the Gromov-Witten invariants of a point. A striking theorem of [19] relates, via 
the method of symplectic reduction, the Weil-Petersson volume of M. g , n {L) and the 
intersection numbers involving both the first Mumford class k\ and the V'-dasses 
on M g , n - As a consequence, she proves that the volume Vol {M. g ^ n (L)'), after an 
appropriate normalization with powers of ir, is a polynomial in L with rational 
coefficients. 

Since there is no particular reason to believe that there should be a direct relation 
between the Weil-Petersson volume of the moduli spaces of bordered Riemann 
surfaces and Virasoro constraint condition, our discovery suggests the existence of 
another, more algebraic, point of view in the Mirzakhani theory. 

Our Virasoro structure also bears an interesting consequence: it leads to the 
natural normalization of the Weil-Petersson volume of the moduli spaces of bordered 
(or unbordered) Riemann surfaces. Although the geometric orbifold picture and the 
algebraic stack picture give the same moduli space for most of the cases, there is one 
exception: the moduli space of one-pointed stable elliptic curves A^i.i- If we define 
this space as an orbifold, then its canonical Weil-Petersson volume is £(2) = n 2 /6. 
On the other hand, the Virasoro constraint condition dictates that we need to have 

Vol(AT M ) = ^ 

as its canonical symplectic volume. This makes sense if we consider Aii.i as an 
algebraic stack. The factor 2 difference is due to the fact that every elliptic curve 
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with one marked point possesses a Z/2Z automorphism. It is remarkable that even 
a purely hyperbolic geometry argument leads us to this stack picture. 

To summarize our main results, let us consider the rational volume of M. g , n {L) 
defined by 

m def Vol (M g , n (L)) 

Vg,n{-Li) - 2 d lT 2d 

where d = 3g — 3 + n and 



Vol(M g , n (L))= f 

JM g , n 



WP 



is the Weil-Petersson volume of At 9in (L). Then the Mirzakhani recursion formula 
reads 

fL\ poo poc 



2 /■ l j 1 poo f- oc 

v g ,n{L) = — / / xyK(x + y,t)v g -i t „ + i(x,y,L i )dxdydt 
J^i Jo Jo Jo 

2 rLi roc r-oo 

+ T E / / / xyK{x + y,t)v guni (x,L x ) 

Ll n.TZ=n JO JO JO 



91+92=5 
II]J={2 n} 

x v g2yn2 {y,Lj)dxdydt 
+ r E / x(K(x,t + L 3 )+K(x,t-L 3 )) 

1 j=2 J° J° 

x Vg tn -i(x,L{Qy)dxdt, 
where the kernel function of the integral transform is given by 

K(x,t) ! ! 



1 + e w ( x+t ) 1 + e*^- 1 ) ' 

and the symbol ~ indicates the complement of the indices. Recall that our nor- 
malized Weil-Petersson volume is a polynomial in L with coefficients given by in- 
tersection numbers of K\ and V'-classes: 

n -. oo 

v*m= e n^n^iw* 

doH hd„ i=0 1 i=l 

=d 

Instead of defining our generating function directly from these rational volumes, let 
us consider the generating function of the mixed k\ and ?/>-class intersections 

a 

= E E («v-^n^- 

9 m,{rii} i=0 

The main results of the present paper are the following differential version of the 
integral recursion formula. 

Theorem 1.1. For every k> —1, let us define 

lA, n , . . . , d l^(2(j + fc) + l)!! d 
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di+d2=k— 1 
di,d 2 >0 

where a, L = ^i+iy. ■ Then we have: 

(1) TVie operators Vk satisfy Virasoro relations 

[V„, V m ] = (n - m)V n+m . 

(2) The function exp(G) satisfies the Virasoro constraint condition 

VfeCxp(G)=0 fork>-l. 

Moreover, these properties uniquely determine G and enable one to calculate all co- 
efficients of the expansion. Since G contains all information of the rational volumes 
v g _ n {L), we conclude that the Virasoro constraint condition is indeed equivalent to 
the Mirzakhani recursion relation. 

Since the generating function for V'-class intersections 

n*d,*i,-) = E< eSTit4 >. 

9 

= £E<IK%nS 

9 {«,} 

is a solution of the KdV hierarchy, it is natural to ask if G satisfies any integrable 
equations. Indeed, we prove the following. 

Theorem 1.2. The function exp(G) is a r-function for the KdV hierarchy for any 
fixed value of s. In fact, we have an explicit relation 

(1.1) G(s,t ,ti,...) = F{t ,t 1 ,t 2 +72,^3 +73,- ••). 

where 7* = (fc^ji!^" 1 - 

We remark that it is well known to algebraic geometers that generating functions 
F and G contain the same information [2, 3, 7, 12, 16, 26]. 

An important consequence of Theorem 1.2 is that G is also completely deter- 
mined by the property of being a r-function, together with the string equation 
V-i exp(G) = 0. It is fruitful to think of the string equation as being the initial 
condition for the KdV flow. Since G is determined, we note that Theorem 1.2 is 
again equivalent to Mirzakhani's recursion formula. 

Here we recall that in the theory of integrable systems, every variable has a 
weighted degree so that all natural operators have homogenous weights. Coming 
from the KdV equations, we assign deg tj = 2j + 1 . The quantity jj has the same 
degree, which defines that degs- 7 = 2j + 3. The Virasoro operator Vk then has 
homogenous degree — 2k for every k > — 1. Another way to view the degree of s l 
comes from the generalized Kontsevich integral 

log / e ^i=° K 2> 3 e 2 dX. 

JHn 

As indicated in the work of Mondello [20] , there should be a substitution Sj = cjs^ 1 
which transforms the asymptotic expansion of the integral into the generating func- 
tion G. Since Sj has degree 2j + 1, we confirm that s- 7 must have degree 2j + 3. As 
well, we should point out that it is quite natural for (1.1) to leave variables to and 
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t\ unchanged. The reason is that in any expression relating intersections involving 
k classes to those involving r terms alone, To and n never make an appearance. 
A few of the natural questions that crop up from this work are: 

(1) Is there a matrix integral expression for the function G? 

(2) Is it possible to prove that G is a solution to the KdV hierarchy without 
appealing to the Witten-Kontsevich theorem? 

(3) What is the direct geometric connection between the cut and join mecha- 
nism and the Mirzakhani recursion? 

We note that the essence of the original Virasoro constraint conjecture is that 
the generating function of Gromov-Witten invariants should have a matrix integral 
expression. The analysis of matrix integrals [1] indicates that once a matrix integral 
formula is established, the Virasoro constraints and integrable systems of KdV type 
are obvious consequences. Since the ribbon graph expansion method provides a 
powerful tool to matrix integrals, the very existence of both the KdV equations 
and the Virasoro constraint for the Weil-Petersson volume of the moduli spaces of 
bordered Riemann surfaces points to a matrix model expression and ribbon graph 
interpretation of the Mirzakhani formulas. These questions, however, are beyond 
the scope of our present work. 

This paper is organized as follows. In section 2 we review the work of Mirza- 
khani [18, 19]. Since the Virasoro structure very delicately depends on all the subtle 
points of the theory, we provide a detailed discussion on some of the key ingredients 
of the work, including the case of genus one with one boundary, precise combinato- 
rial description of cutting a surface along geodesies, and the choice of a canonical 
orientation of the circle bundle when the Duistcrmaat-Hcckman formula is applied 
to the extended moduli spaces. Section 3 gives a proof of Theorem 1.1. Finally, in 
section 4 we prove Theorem 1.2. 

The second author would like to thank Greg Kuperberg and Albert Schwarz for 
helpful conversations. 

2. Mirzakhani ' s Recursion Relation 

2.1. Notations. Since the orbifold picture and the stack structure are the same 
for the moduli spaces of algebraic curves except for genus 1 with one marked point, 
we employ the orbifold view point throughout the paper. As mentioned above, 
however, when we interpret the canonical volume of the moduli spaces, we need to 
use the stack picture. 

Let M.g y n denote the moduli space of smooth algebraic curves, or equivalently, 
the moduli orbifold consisting of finite area hyperbolic metrics on a surface, of 
topological type (g,n). A surface of type (g,n) is a surface with genus g and 
n punctures. Since we are interested in the stable, noncompact case, we impose 
2g — 2 + n > and n > throughout this paper. When referring to the underlying 
topological type of the surface we will consistently employ the notation S gin . We 
will also use the notation Ms for the moduli space of surfaces of topological type 
S. 

Mirzakhani's breathtaking theory is about the moduli space M. g ^ n (L) of genus 
g hyperbolic surfaces with n geodesic boundary components of specified length 
L = (Li, . . . ,L n ). This space relates to the algebro-geometric moduli space via 
the equality M. g , n — -M s>n (0). The moduli space of bordered Riemann surfaces is 
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defined as an orbifold 

M g , n (L) = T Si „(L)/Mod Si „, 

where Mod 9j „ is the mapping class group of the surface of type (g, n), i.e., the set 
of isotopy classes of diffcomorphisms which preserve the boundaries setwise, and 
Tg :n (L) is the Teichmiiller space. The Deligne-Munford type compactification of 
this moduli space is obtained by pinching non-trivial cycles. 

The tautological classes we consider in this paper are the m and tp classes. Let 

7T : Mg, n +l{ = Cg tn ) ► Mg, n 

be the forgetful morphism which forgets the n + 1-st marked point, and 

(Ji{C,xi,. . .,x n ) = XiE C, i = 1,2, . . . ,n 

its canonical sections. We denote by luc/m the relative dualizing sheaf, and let 
T>i = a; L {M. g n ), which is a divisor in M g n +i- The tautological classes are defined 

by 

Ci = ct*(ujc/m), 
tpi = ci(A), 

«i = 7T* ^ci(w c/-m (x;c(a)) 2 

We arc interested in the intersection numbers 



«r dl ---r dn ) g = [_ 

JMg,„ 



The class K\ has a nice geometric interpretation, coming from the symplectic 
structure of M. g , n - The Fenchel-Nielsen coordinates are associated with a pair of 
pants decompostion of the surface S g>n , which is a disjoint set of simple closed 
curves T = {71, . . . ,7^} such that S g , n \ T is a disjoint union of pairs of pants 
(triply punctured spheres). Since pairs of pants have no moduli (they are uniquely 
fixed after specifying the boundary lengths) , all that remains to recover the original 
hyperbolic structure is to specify how the matching geodesies are glued together. 
Hence for every curve 7$ in the pair of pants decomposition, one has the freedom 
of two parameters (U,Ti) where U is the length of the curve and r» is the twist 
parameter. These coordinates give an isomorphism 7^„ = x R d . At the moduli 
space level we must quotient out by the mapping class group action. Since one full 
twist around a curve is a Dehn twist (element of Mod 9! „), we get local coordinates 
of the form R d + x (S ,1 ) d . 

In Fenchel-Nielsen coordinates, the Weil-Petersson form is in Darboux coordi- 
nates [24]: 

uy/P = / dli A d,Ti. 



This is a closed, nondegenerate 2-form on T g ^ n which is invariant under the action 
of the mapping class group, hence gives a well defined symplectic form on -M ffi „. 
Note that by Wolpert [25] the Weil-Petersson form extends as a closed current on 
M.g, n - In particular, the Weil-Petersson volume 



Ja 



WP 



M g , n (L) d\ 
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is a finite quantity. (Note that Wolpert defines the Weil-Petersson form to be half 
of the above expression. Our convention is adopted from the algebraic geometry 
community [2, 12, 26].) The relation to tautological classes is provided by the 
well-known formula 

u>wp = 27T K\. 

For use in the sequel, we note that the vector field generated by a Fenchel- 
Nielsen twist about a simple closed geodesic is symplectically dual to the length of 
the geodesic. A Fenchel-Nielsen twist is defined by cutting the surface along the 
curve, twisting one component with respect to the other and than regluing. As a 
formula, we have 

wwp( ■ , tt- ) = dk. 
on 

2.2. McShane's identity. A crucial step in Mirzakhani's program [18, 19] is to use 
McShane's identity to write a constant function on the moduli space as a sum over 
mapping class group orbits of simple closed curves. To state Mirzakhani's general- 
ization of McShane's identity, we introduce the following notation for an arbitrary 
hyperbolic surface X with boundaries {j3\, . . . , (3 n ) of length (L\, . . . , L n ): Xj de- 
notes the set of simple closed geodesies 7 such that (f3\ , (3j , 7) bound a pair of pants; 
and J the set of pairs of simple closed geodescis (0:1,0:2) such that cti, o 2 ) 
bound a pair of pants. Using the functions 

/ e x/2 + e (y+z)/2 \ 

T>(x,y,z) = 2log I /2 + e(y+z)/ \ and 



/ cosh I + cosh 
K cosh f + cosh ^ 



Mirzakhani proves 

Theorem 2.1 (Mirzakhani [18]). For X as above, we have 

n 

Li= -D{LiA^i)A^))+Y.H n ^ L ^ 1 ^)- 

2.3. Integration over the moduli spaces. The idea is to find a fundamental 
domain for a particular cover of M g , n , enabling one to integrate functions over this 
covering space. Then a specific class of functions defined on the moduli space (such 
as those arising from McShane's identity) can be lifted to this cover and integrated. 
To that end, let T = {-fi, . . . , 7„} be a collection of disjoint simple closed curves on 
the surface S g , n , where S g , n is the underlying topology of the hyperbolic surface 
X e Mg tn . We define 

StabT = n Stab 7 , = {/ G Mod s ,„ | = 7i , for i = 1, . . . , n}, 

and set 

-^L = r ^«/ Stabr 

= {(A, 771, . . . , r] n ) I X e M- g ,n, Vi is a simple closed geodesic in Mod7i }. 

Note that as a quotient of Teichmiiller space, A4^ n inherits the Weil-Petersson 
symplectic form. Hence we can talk about integration over n with respect to 
the symplectic volume form. The advantage of integration on A4g „ as opposed to 
the usual moduli space is that we can exploit the existence of a hamiltonian torus 
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Figure 1. Decomposing a surface 



action. In fact, by a result of Wolpert, the vector field generated by a Fenchel- 
Nielsen twist along a geodesic is symplectically dual to the length function of the 
geodesic (as a function on Teichmiiller space). The space M r g n is the intermediate 
covering space on which the circle actions on {71, • • • , 7„} descend. The problem 
with attempting to construct such a circle action on M g , n is that there is no well 
defined notion of a geodesic curve on an clement of moduli space; the best one can 
obtain is a mapping class group orbit of curves. 

Thus we have the moment map for the torus action: 

I : M T g>n - R™ 
(X,r,)~ (l( Vl ),...,l(r] n )). 

Hence we see that l~ 1 (x)/T n is symplcctomorphic to Ms g n \r(L, x, x). Recall 
that Ms gn \r is the moduli space with underlying topological type the (possibly 
disconnected) surface 5 Si „ \ T and with boundary lengths following the rule outlined 
in Figure 1. 

The most straightforward way to prove the above assertion is to take a pair of 
pants decomposition for the surface S gt „ which contains the curves T. Then/ -1 ^) 
fixes the lengths of the geodesies 771 , . . . , t] n , while quotienting by the torus action 
removes the twist variable from these curves. This gives the diffeomorphism. The 
symplectic equivalence follows immediately from the Fenchel-Nielsen coordinate 
expression for the Weil-Petersson form. What emerges is an exceptionally clear 
local picture for the space M T g n . In fact, it is a fibre bundle over R™ where the 
fibres are (locally) equal to the product of a torus and Ms g n \r(L, x, x). 

Consider a map 

/ : M T gtn R, 

which is a function of the lengths of the marked geodesies l(rii). In other words, 
f(X,r]i, . . . ,rj n ) — /(Z(?7i), . . . ,l(r] n )). By the previously discussed decomposition 
of A4^„ we can write 

(2.1) f f(l(v))e" WP ^ = [ f( x )Vol Sgn \ T (L,x,x)x-dx. 

JM T gn JRl 

Here e ulWP means we are integrating over the maximal power of the Weil-Petersson 

uj d 

form — , d = 3g — 3 + n. Note that if 5 9i „ \ T is disconnected then M.s g „\r is the 
direct product of the component moduli spaces, with the volume being the product 
of each. 
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2.4. Volume calculation. To relate the above discussion to integration on the 
moduli space, Mirzakhani uses her generalized McShane identity. As a lead in 
to the main results, consider the following simplified situation. Suppose 7 is a 
simple closed curve on S g , n , with Mod 7 its mapping class group orbit. Then given 
any hyperbolic structure X on S g , n , every a G Mod 7 has a unique geodesic in 
its isotopy class. Denote lx{ot) the corresponding geodesic length. Hence for any 
function / : R + — > R and thinking ofle [X] as a representative of an element of 
A4g in we have the following well defined function on A4 gin : 

P ■ M g , n - R 

a£Mod 7 

One can easily check that this function docs not depend on the choice of repre- 
sentative X g [X]. However, it is not a priori clear whether or not p will be a 
convergent sum. At minimum one requires lim^^oo f(x) = 0. We can similarly 
define a function 

f ■■ Ml n - R 

by the rule 

P(X,r,) = f(l x (r,)), 

which gives the relation 

p(x)= ]T P(x, v ). 

In particular, since the pullback of the Wcil-Pctcrsson form is the Wcil-Petersson 
form on the cover, we have 

j-J e ^wp _ / jl e ^wp 

Mg,n JMl, n 

Note that for any curve 7 g Ij, we have Tj = Mod 7. This follows because two 
curves are in the same orbit of the mapping class group if and only if the surfaces 
obtained by cutting along the curves are homeomorphic, with a homeomorphism 
preserving the boundary components setwise. The homeomorphism will extend 
continuously to the curves to give the map of the entire surface. This tells us that 
the set J is not the orbit of a single pair of curves (ai, 012) S J ■ In fact, we further 
refine this set of curves as follows. For any (ai, 0:2) € J set P((3i, a\, 0:2) C 5 fl ,„ to 
be the pair of pants bounded by the curves /3i, ai, 0:2- Now we define (see Figure 2) 

Jconn = {(ui,a 2 ) e J I Sg tn \ P((3 1 ,a 1 ,a 2 ) is connected } 

»7gi,{ii,...,* ni } = {("1^2) g J I Sg tTl \ P(f3i,a 1 ,a 2 ) breaks into 2 pieces 

one of which is a surface of type (gi, ri\ + 1) with 

boundary (a*, /3 ix , . . . , /3 ini ) }. 

Other than the obvious identification 

<Jgi,{ii,-..,i ni } ~~ "9—91, "}\{»iv>»ni}' 

these subsets form disjoint orbits under the mapping class group. Moreover Mod giU 
acts transitively on each set. 
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Figure 2. Removing a pair of pants from a surface 



Hence we can write the Mirzakhani-McShane identity in the following form 

L ^ = \ E E V(L u l( ai ),l(a 2 )) 

31+32=3 (ai,Q 2 )eJ sl ,A 

{2,...,n} 

+ Yl 'D(L 1 ,l(5 1 ),l(5 2 )) 

(5l,52)£Jcon„ 

n 

+ EE^( L i<^(7))- 

j=2 7 ei 3 - 

There is a slight inaccuracy - we undercount by half for terms with n = 1 and 
.9i = .92- However, we will see in a moment that this makes further calculations 
somewhat simpler. 

Each of the terms in the above sum can be lifted to a function on an appropriate 
cover M^„. We see that 



/ l^ p = \ EE/ 

JM g , n (L) * a . +a , = a( n , _ n ^<=v „ JA 



91+32=3 ( ai ,a 2 )eJ gi ,A JM 9.^) 
A\\B= 
{2,...,™} 



V{L u l{ ai ),l{a2))e u 



Hence 



+ Y, / ©(Li,/(5 1 ),Z(<5 2 ))e^ 

(5i,5 2 )GJco»„ JM ^n{L) 
n „ 

+ EE / nL u L h l{l))e" WP ■ 

3=2 1 £l ] J M g , n (L) 

L 1 \o\ g , n {L) = \ J2 [, ML u l{ m ),l{m)) 
2 3i+^=3 <>MlT a2} 
A]\B={2,...,n} 



2 JMl'i-'** 



v(L u i(m)Am))^ wp 

Y / ^(L 1 ,L J -,/( 7 ))e^. 
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Note that the factor \ that appears in front of the second term in the sum 
is needed because the mapping class group orbit double counts the set of curves 
(61,62)- In other words, there is a diffcomorphism exchanging 61 and 62- Similarly, 
the factor i discrepancy for the term in the first sum with gi = g 2 and n = 1 has 
now disappeared and the presented sum is unambiguously correct. Applying the 
results of the previous section we see that 

Li \o\ g , n {L) = ^ Y f xyViL^x^) Vol guni (x, L A ) Vo\ g2 . n2 (y, L B )dxdy 
ai+g2=g jR + 

+ \ f xyV(L 1 ,x,y)Vo\ g ^ 1 . n+1 (x,y,L i )dxdy 
1 Jm 2 + 

+ ^ / x1l(Li_,Lj, x) Vo\ g ^ n -i(x,LQ)dx. 

j=2 J& + 

For any subset A = {ii, • • • ,lk} C {1, . . . , n} the notation L A means the vector 
(L h , . ..,L lh ) while L A = i{i ; ...,„}\^. 

There is one additional subtlety that crops up at this point. Note that for the 
case of A / Ii i i(i), there is an order two automorphism obtained by rotating around 
the boundary by half a turn. There are two ways to deal with this issue. The 
approach taken in [18] is to divide the appropriate integrals by 2 every time such 
a term appears in the above integral. Our approach, which is computationally 
equivalent, is to define the volume of A4u(L) to be half the value obtained by 
calculations using the above techniques. In other words, we have initial conditions 

Vol , 3 (£) = 1 

Vol 1 , 1 (L) = l(L 2 + 4 7 r 2 ). 

We will see that this viewpoint simplifies further calculations; as well, it agrees with 
known results from algebraic geometry. 

The final step is to differentiate both sides with respect to L\ and then integrate, 
which has the effect of simplifying the integrands on the right side of the equation. 
The result is 

fL-\ r-OO fOC 



1 r^i /"oo poc 

Vo\ g , n (L) = — J2 xyH(t,x + y) 

fi 1 — U — n 



91+32=3 ' 

x Vol Sl >ni (x, L A ) Vol 92 ,„ 2 (y, L B )dxdydt 



OO />OC 



where 



x Vo\ g -^ n+ i(x,y,Li)dxdydt 

+ ^rE / / x ( H ( x > L i + L >) + H ^ L i - L i)) 

x Volg in -i(x,LQ)dxdt, 

H ( X > V) = l +e (x+y)/2 + l +e (x-y)/2- 
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Figure 3. Capping off a bordered surface 



2.5. Relation to intersection numbers. In this subsection we review the idea 
of Mirzakhani to write the integral over A4 g , n (L) as an integral of an appropriately 
modified volume form over M g ^ n - This will relate the Weil-Petersson volumes to 
intersection numbers of tautological classes. Following [19], let 



be the moduli space of bordered hyperbolic surfaces of arbitrary boundary length, 
with the additional information of a marked point on each boundary component. 
If Li = 0, then we can think of pi as a point on a horocycle about the cusp. 

The marked point can be used as a twist parameter, so by gluing on pairs of 
pants with two cusps and the third boundary having length matching the surface's 
boundary, we obtain a map M. g , n — * M.g,2n- I n fact, we have M g , n = M. g 2n where 
r = {71, . . . , 7„} is a collection of curves which group the cusps into pairs. We refer 
to Figure 3 for a descriptive picture of this construction. 

This tells us that M. g ,n has a symplectic structure from the Weil-Petersson form 
on MF 2 n - Moreover it has a hamiltonian torus action given by rotating the marked 
points on the boundary. However, we need to take some care here. We are interested 
in studying the symplectic action in a neighborhood of surfaces X G M g . n with 
— 0. Moreover, we want to construct the action in such a way that these 
points are not fixed by the torus action. In other words, we need to non-trivially 
extend the action to the cusped surfaces. 

It is a simple matter of defining the twists to be proportional to the lengths of 
the boundaries. In other words, we scale the action so that a twist parameter of 1 is 
always the identity. The model is the change from the cartesian (x, y) coordinates 
in the plane to the polar coordinate (r, 6). In the first case rotation around the 
origin leaves it fixed, but (0, 9) is not fixed by 9 9 + e. From the point of view 
of Mg^n, the marking on the boundary degenerates to a marking on a horocycle 
of the cusp. The result, after a change of coordinates to the reparametrized twist 
coordinate 9- L — Ti/li, is 



X e Mg.n to the lengths of its boundary components, we see that L~ (x) is a 
principal torus bundle over M. g>n (x). In fact, over Ai gt n(0) it is the principal 
bundle associated to the vector bundle C At first glance this is a rather 



M 



3," — 



{(X, P1 , . . . , Pn ) I X e M g , n {L), L e R| ,p* e ft} 




MIRZAKHANFS RECURSION AND VIRASORO CONSTRAINTS 



13 



counter- intuitive statement, as marked points on the boundary map naturally to 
the tangent bundle, rather than the cotangent bundle. However, the principal 
torus bundle in question is naturally oriented from the induced orientations on the 
boundaries coming from the orientation of the surface. This orientation is opposite 
to the natural complex orientation on the tangent bundle. This is most easily seen 
by studying the clockwise orientation induced on the unit circle from the standard 
orientation of the plane. 

Using symplectic reduction, we see that the reduced space L~ 1 (x)/T n is sym- 
plectomorphic to A4 g , n (x) with the Weil-Petersson form. We may use the tech- 
niques of the Duistermaat-Heckman theorem to compare u>wp(L) to u>wp(0). The 
result is 

u) WP (L) = u WP (0) - - L i Curv(A), 
where Curv(£j) is the curvature of the bundle. Since ci(£j) = — Curv(£i) we get 



lj W p(L) = u WP (0) + i E Lfyi- 



2.6. A rational recursion relation. Using Wolpert's equivalence m — we 
define the rational volume of M. g , n {L) to be 

(T . def Vol 9i „(27rX) 
v aA L ) = h-^r d d = 3g - 3 + n 



2 d TT 2d 

d\ 



n 1 oo 

= e n^n^un^- 



<k- 

doH hd„ i=0 »=1 

=d 

We reformulate Mirzakhani's recursion relation for Vol ff , n into a recursion relation 
for v g , n . Making the above change of variables to the recursion relation gives 

2 f'^i r°° r°° 

v g , n (L) = — / / xyK(x + y,i)v g - lin+ i(x,y,L{)dxdydt 

Jo Jo Jo 



2 

i 



+ — V / / / xyK(x + y,t) 

Ll 91+92=9 J ° J ° J ° 



OO POO 



x v gi , ni {x, Lx)v g2 >n2 (y, Lj)dxdydt 
+ T E / / x ( R ( x > 1 + L i) + K ^ 1 - L ^ 

1 j=2 ^° J{) 

x v gi „-i(x, L^Qy)dxdt, 

with normalizations v ^(L) = 1 and v\.\{L) = ^(1 + L 2 ). 
The integral kernel K{x,t) is defined as 

K i X ^) = I + e7r(x+t ) + 1 + g W (x-t) ' 

which gives the following integral identities: 

a t r°° r 2k+1 

(2-2) fok+iit)^^ -^—yK(x,t)dx 
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fe+1 



p> j-2fe+2-2m 

E/_1 \m-l(r,2m _ q\ D 2m t 
_ n l ' { >(2m)\ (2fc + 2-2m)!' 



m=0 

/•oo /-oo x 2i+l 2j + l 
(2.3) f) 2i+2j+3 (t)=y o ^ _____^( X + y,^y. 

Note that -B 2m is the 2m-th Bernoulli number. 
3. From Mirzakhani's Recursion Relation to the Virasoro Algebra 



Our aim is to show that the Mirzakhani recursion relations are equivalent to an 
algebraic constraint on the generating function for m and ip class intersections. We 
introduce the formal generating function for all K\ and ip class intersections 

dcf ■ 

a 

3 m,{n;} i=0 

The main result of the paper is the following. 

Theorem 3.1. There exist a sequence of differential operators V-i, Vq, V\, . . . sat- 
isfying Virasoro relations 

[V n ,V m ] = (n- m)V n+m 

and annihilating exp(G): 

Vfeexp(G)=0 for k = -1,0,1,... 

This property uniquely fixes G and enables one to calculate all coefficients of the 
expansion. 

The proof is obtained by differentiating Mirzakhani's recursion relation. For 
reference, we note that 

^2fc! £)2fc„ 

<3 - 1) m^--aSF v ' AL) 



di>ki 



and 



(3 - 2 » ^•••^•Bn(f)<"N-^ 

where fco = Zg — 3 + n — 5TJ™=i 

The recursion relation gives the following identity for (g, n) ^ (0, 3), (1,1). 



2fe„ 



d 2fcl & 

Q2ki 2 f Ll f°° f°° 



dL 2k 



fei 2 f 1 r°° r°° 

21777/ / / x v K ( x + y^) 

! L l JO JO JO 
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Q2k 2 Q2k n 

x dLf 2 ' ' ' dLf- 

d 2kl 2 ^ f Ll f°° f°° Tjr . 



v g -i,n+i(x,y, L^dxdydt 



L=0 



gi+92=g 



Q2k(I) Q2k(J) 

x 2k{T) v 9i^i ( x ' L z) 2~kjj) v 92,n 2 {x, Lj)dxdydt 



L=0 



" Q2(k 1 + k ] ) i too 



dL 2k(l,j) 
1,3 



v g ,n-i{x, LQ)dxdt 



L=0 



Plugging in the expressions for derivatives of volume functions (3.1), (3.2) and 
integrating against the kernel function using (2.2) and (2.3) gives 



v (2rf 1 + l)!(2rf 2 + l)! A (2fcQ! d 



d + d l + d 2 

= fc +fc 1 -2 



71+ 1 



a 2A;i 2 r^ 1 

^TJ ^(*-Hfa)-H.(*)* 



Li=0 



v v (2rf 1 + l)!(2rf' 1 + l)! A(2fej)l , d 

+.92=fld +di=3oi-3 u i o 1 l=2 



'31,11 



gi+92=g d +di=3gi-3 
XIJJ +m-fc(I) 

d +d 1 =3g 2 -3 
+n 2 -k(J) 



E E 

j=2 d +di = 
feo + fci+fej — 1 



d' <9 2A:i 2 f^ 1 

rfo!rfl! 11 fcl \ K 1 r * T fc(lJ)/ fl ,n-l 



Li=0 



a 2fci 2 /" Li 



5Lf 1 ii 7o 



Li=0 



We rewrite this sum by introducing the sequence of nonnegative integers {n , 
m, «2, • ■ •} such that rij = # \{ki \ i ^ 1, k% = j}\, and relabel k\ = k. In other 
words, we have 

(^V fel ...r fen ) s = (^V fe r "Vr---) s . 
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We further define 

(3-3) ^ = (-1)^(2^-2)^, 

which results in the equation 



(3.4) (2fc + l)!!(^r t []rr) 9 

i=0 

1 - h I °° 

= 2 E fi%o^o)(2di + 1)1(2* + l)!(^r dl r d2 nrD s -i,„ +1 

d +di+d 2 = i=0 
feo+fc-2 

+ ^ E E #T%o- d o-4)^ + 1 ) ! ( 2d i + 1 ) ! 

31+32=3 d +di=3gi-3 u ' °' 
{i,} + {m,} = {n,} +ni-fe(X) 

d +d 1 =3g2-3 
+n 2 -fe(J) 

oo i , 



i=0 



+E E ^-J^^M^r^).- 

j=Q d„+d 1= u ' v J >' 

ko+k+j-1 



By looking at expressions of the form 



5 m^rii} i=0 



at 



^.g -EE ^<»rvw - ■ •>.£ ft £. 

3 m,{ni} 'i=0 

d 2 G ml s m °° 

Sl 9^ = E E (^)!^ r ~ VjTfeTo " OTl " 1 "'^11^!' 

9G 9G ■» ■» ml s m 00 t™' 

9 31+32=3 Z i=0 

"M";} di+d 2 =m-i 

{ki} + {h} = M 

we see that (3.4) leads to the following expression for all k > 0: 

i+fe-l 

Note that similar expressions are possible for k = — 1, by taking special care of the 
base cases (g,n) — (0,3), (1, 1). We introduce the family of differential operators 
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for k > — 1 

2 9t fe+1 ' u 4 48^ 48 

1 ~ {2 (i+j + k) + iy.\ a 



2 ^o ( 2 ^'- 1 ) !! J dU +J+k 



3 ^H+j+k 

d 2 



«— di+d2 — 
i+fe-1 

We have proven the following statement. 
Theorem 3.2. For k > -1 

14 cxp(G) = 0. 

A reasonable question is: what is the algebra spanned by the operators V^l The 
answer is that they span a subalgebra of the Virasoro algebra. One can check 
directly that the operators satisfy the relations 

oc 

[V n , V m ] = (n-m)^2 /3 4 s 4 K +m+4 . 

i=0 

On the surface, this looks to be a deformation of the Virasoro relations (setting 
,s = recovers Virasoro). However, this is, in fact, a simple reparametrization 
of the representation. These statements can all be proved by direct calculations. 
Here we make some simplifications. Let us introduce new variables {T 2 j+i}j=o,i,... 
defined by 

U 

T2l+1 = (2i + l)!!' 
which transform the operators Vk into 

k 2dT 2k +3 • U 4 48 ; 16 



+ \ E ( 2 J + W^+i 

i,j=0 
oo 

+ iE E a 



2 ~f ' d T 2(i+j+k)+i 



i+k-l 

This admits the following 'boson' representation, similar to that used by Kac and 
Schwarz [11]. Define operators J p for p e Z by 



(-p)T_„ ifp<0, 



" ^ a if p >o. 



Then 

1 

V k = --J 2 fe+3 + ^fos'Ek+i, 



i=0 
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where 



Ek — J E ^2p+1^2(fc-p)-l + 



16 ' 



To recover operators satisfying the Virasoro constraint we need a better handle on 
the constants as defined by (3.3). Starting from the defining formula for the 
Bernoulli numbers 

z/2 



B 2n 2n _ £ £^_+e_ 

^ (2n)\ Z ~ 2 e*/ 2 - e-*/ 2 ' 



we see that 



2s(cot V«/2-cot V2s) 
2s 



»=0 



sin V 2s 

This motivates the definition of the constants by the series 

sin \[2s 



E 

i=0 



QfiS 



2s 



from which we obtain the operators 

oo 

Vk = E a * st ^fc+* 

*=0 
^ oo 

(3.5) =-- VaisVsfe+s + ^fe. 

We are now ready to prove the following. 

Proposition 3.3. The operators Vk, k > — 1 satisfy the Virasoro relations 

[V„, V m ] = (n - m)V n+m . 

Proof. The first step is to verify that operators Ek satisfy the Virasoro relations, 
which is a straightforward calculation. Since \J2k+3,E m ] = (2k + 3)J2(fe+ m )+3 we 
see that 



[V„,V m ] = 



1 ' x | 

2 E a i S ^2(n+j)+3 + E n , —- QiS* J 2 ( m+ j)+3 + -E„ 
"2E a ' S '(^ 2 ("+')+ 3,Em ] + [En, J2(m+i)+3]^J 



i=0 

(n - m)V n+m . 



a 
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4. Relationship to KdV hierarchy 

The Witten-Kontsevich theorem [23, 15] states that the generating function for 
tp class intersections 



F(io,*i,...)=E< eI>iti > fi 



=EE(nnn£ 

9 {n,} 

is a T-function for the KdV hierarchy. The property of being a r function, combined 
with the string equation 



<-on^> s =£<n^-> s > 



' 9 

i—1 j — 1 i—1 

completely determines the function F. Another way of determining F is the Vira- 
soro constraint condition. Let us define the sequence of operators Lk for k > — 1: 

UU t (2fc + 3)!! d , 1 " (2(j + fc) + l)ll d 
[ ' k ~ 2 dtk+1 + 2^ (2j-l)!! * 3 'at j+fc 

4, J, Wl + l)H(«. + l)H s ^ + ^ + ^. 

di+d 2 =k — 1 
di,d 2 >0 

The Witten-Kontesevich theorem, together with the string equation, implies 

L fe (expF) =0 

for k > —1. This property is also sufficient to uniquely fix F. Note that L_ie F = 
is equivalent to the string equation. The consistency of the infinite set of differential 
equations follows from the fact that operators L n satisfy the Virasoro relations: 

[L n , L m ] = (n - m)L n+m . 

Recall the operators V k defined in equation 3.5 (rewritten in terms of the variables 

U) 

1^,,„, „.,, , d 1^ (2(j + fc) + l)!! d 

- "5 1> + l> + ^ + 5 § W ^ 

+2 



d 1 +d 2 =k-l ai d2 

d u d 2 >0 

where a» = ^T+i)! • ^he change of variables 



U = 



ti for i = 0, 1 , 

— (2? — l)!!a i _is J_1 otherwise, 

transforms the operators Vk into 



Vfe - " 2 (2 * + 3 )'fe + 2 ^ n (2,-1)!! ^ 



20 



MOTOHICO MULASE AND BRAD SAFNUK 



di,d 2 >0 

But these are precisely the operators L k (4.1). We have thus proven the following. 
Theorem 4.1. 

G(s,to,t lt ...) = F(t ,t 1 ,t 2 +nf2,t3 +73, • • ■), 

where ji — j^jj^s 1 _1 . In particular, for any fixed value of s, G is a r function 
for the KdV hierarchy. 

That the more general generating function G is expressible in terms of F is not a 
surprise. It has been known since at least the work of Witten [23] that intersections 
involving k classes are expressible in terms of i\> classes. Moreover, Faber's formula 
[7] for this correspondence gives an explicit proof of the above theorem. In fact, 
one has 

K n - k -r. . C/.I71I+1 . . ,| 7fc |+ix 

aes n Ili=i(l7i|-1)! 

(o-=7i"'7fe 
is cycle dccomp) 

which gives a short, direct proof of Theorem 4.1. This is essentially the approach 
taken by Zograf [26] for his calculation of the Weil-Petersson volumes of A4 g , n . 
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